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In order to study the S-nucleus optical potential, we measured inclusive {-k~,K'^) spectra on 
medium-to-heavy nuclear targets: CH2, Si, Ni, In and Bi. The CH2 target was used to calibrate 
the excitation energy scale by using the elementary process p -\- tt" K'^ -I- E~, where the C 
spectrum was also extracted. The calibration was done with ±0.1 MeV precision. The angular 
distribution of the elementary cross section was measured, and agreed well with the previous 
bubble chamber data, but with better statistics, and the magnitudes of the cross sections of the 
measured inclusive {■k~,K^) spectra were also well calibrated. All of the inclusive spectra were 
found to be similar in shape at a region near to the S~ binding energy threshold, showing a weak 
mass-number dependence on the magnitude of the cross section. The measured spectra were com- 
pared with a theoretical calculation performed within the framework of the Distorted Wave Impulse 
Approximation (DWIA). It has been demonstrated that a strongly repulsive E- nucleus potential 
with a non-zero size of the imaginary part is required to reproduce the shape of the measured spectra. 
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I. INTRODUCTION 

A hypernucleus is the result of implanting of a hy- 
peron(s), such as A, S, S or AA, in an ordinary nucleus. 
Spectroscopic studies of hypernuclei can provide infor- 
mation on the hyperon behavior in a nucleus. One can 
learn about the hyperon-nucleon interaction, and even 
the baryon-baryon interaction within the framework of 
S'[/i?(3), from the hypernuclear structure through the ef- 
fective hyperon-nucleon interaction in a nucleus. So far, 
the spectroscopic studies of A-hypernuclei are in a rel- 
atively advanced stage. The A major shell structures, 
even up to ^^Y or ^°®Pb, have been clearly observed 
with high-resolution spectroscopy in the (Tr+jif"*") reac- 
tion 11,0, 3 , which have revealed the A-single particle 
nature. Namely, the binding energies of A can be well- 
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riCTjroduced by a Woods-Saxon-type one-body potential 
U. In addition, high-resolution Ge detectors together 
with a magnetic spectrometer have successfully observed 
7 transitions involving the A-hypernuclear states split by 
AA'^ spin-dependent interactions, such as the spin-orbit 
and spin-spin forces in several light A-hypernuclei 0, 0| . 

On the other hand, knowledge about and S- 

nucleus interactions is still primitive, although the his- 
tory of S-hypernuclear studies is almost comparable to 
that of A-hypernuclei. Since a strong conversion pro- 
cess of E {YiN — > AA^) takes place inside a nucleus, the 
S-hypernuclear state is expected to be too broad to ob- 
serve. The first claim of the observation of a narrow 
S-hypernuclear state by the Saclay-Heidelberg collabora- 
tion |3] thus generated much experimental and theoreti- 
cal excitement, but was excluded in the later experiment 
So far, many other experiments have been carried 
out on different light nuclear targets using the {K~ 
reactions 0, ^3 • The only E-hypernuclear bound state 
has been established in |,He produced by an in-flight 
{K~ ,'K~) reaction on '*He 0, which was first claimed 
at KEK The observation of this bound state was 

predicted with the effect of a strongly isospin-dependent 
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E-nucleus potential (Lane term) [l^. The existence of 
the Lane term was also suggested in light nuclear sys- 
tems from the systematic difference between the inclusive 
(if-,7r±) spectra T^. All of the other experiments could 
not find any bound S-hypcrnuclear state, and failed to 
give any significant constraint on the E-nucleus poten- 
tial. No data are available in heavier nuclei of A>16. 
One may construct the E-nucleus folding potential based 
on the two-body hyperon-nucleon potentials. Unfortu- 
nately, the hyperon-nucleon potentials based on the one- 
boson exchange model allow several sets of parameters, 
since hyperon-nucleon scattering data are very limited in 
accuracy as well as in statistical quality, as compared to 
the very precise data on the NN system. The informa- 
tion from the E~-atomic data also leaves various shapes 
of possible potentials inside a nucleus, because the X- 
ray data are mainly sensitive to the potential outside the 
nucleus. Information concerning the E-nucleus poten- 
tial can provide a strong feedback not only in nuclear 
physics, but also in astrophysics. The role of strangeness 
in highly dense matter in a neutron star core is being 
intensively discussed 0, Q, 0|. The E~ interaction 
in nuclear matter is of particular importance, since E~ is 
expected to appear first in dense neutron matter to mod- 
erate the chemical potential while maintaining the charge 
neutrality due to its negative charge |l5l |. Therefore, new 
experimental data on medium heavy nuclei are desired in 
order to extract a conclusive understanding about the E- 
nucleus potential. Particularly, the central part of the 
E-nucleus potential is expected to be significant in heavy 
nuclei, where the Lane term is less dominant due to its 
dependence. 



II. PRESENT EXPERIMENT 

The present experiment (KEK-PS-E438) was carried 
out at the K6 beam line of the KEK 12-GeV proton 
synchrotron (PS) We measured inclusive {n^ ,K^) 
spectra for the first time on medium-to-heavy nuclear tar- 
gets: CH2, Si, Ni, In and Bi. Since the inclusive (7r~,_ft'+) 
spectrum reflects the final state interaction of a E^ with 
the residual nucleus, an analysis of the spectral shape 
would give the E-nucleus optical potential. The elemen- 
tary reaction of 

TT- + p K+ + j:- (1) 

was measured with a CII2 target in order to calibrate 
the energy scale and the cross section of the inclusive 
spectra. It is known that the elementary cross section at 
the forward-scattering angle decreases rather smoothly 
with an increase of the incident-beam momentum, while 
the momentum transfer (Ap) increases rapidly below 1.2 
GeV/c 01 . A rapid change of Ap makes the spectrum 
analysis complicated and requires a wider momentum 
acceptance spectrometer system for the scattered kaon. 
The beam momentum was thus chosen to be 1.2 GeV/c. 




FIG. 1: Schematic view of the experimental setup. 



III. EXPERIMENTAL SETUP 

A schematic view of the whole experimental setup is 
shown in figure Q It is composed of two parts: a beam 
spectrometer, used to measure the incident pion momen- 
tum, and the SKS system for the scattered kaon momen- 
tum. A detailed description of the spectrometer system 
can be found in Ref. llSl. 



A. Beam spectrometer 

The beam-line spectrometer comprises a dipole mag- 
net and four quadruple magnets (QQDQQ) together with 
four sets of high-rate drift chambers (BDCl-4), a freon- 
gas Cerenkov counter (eCC) and two sets of segmented 
plastic scintillation counters (BHl and BH2). The elec- 
tron contamination in the beam was rejected by eGC 
with a rejection efficiency better than 99.9%. The BHl 
was segmented into seven vertical pieces of 5 mm thick 
and was installed downstream of eGC, while the BH2 was 
segmented into 5 vertical pieces of 3 mm thick. In order 
to reduce the energy-loss straggling in this counter, it 
was made as thin as possible and placed 40 cm upstream 
from the experimental target. The BH2 was used as a 
time-zero counter for timing measurements of incident 
and scattered particles. By requiring timing coincidence 
between BHl and BH2, any contamination in the pion 
beam was rejected, where the beam trigger is defined as 
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BEAM= BHlxBH2xeGC. 

The beam-line drift chambers (BDCl-4) were placed 
upstream and downstream of the QQDQQ system. In 
order to operate under a high counting rate of several 
M/spill, the sense- wire spacing were made to be short (5 
mm), where the drift space was ±2.5 mm. Each cham- 
ber had six layers of sense-wire plane (xx'uu'vv'), where 
the vertical and ±15° tilted wire planes were denoted by 
x,u and v, respectively. The beam track was measured 
by BDC's with a position resolution of 300 fim in rms. 
The beam momentum was obtained particle-by-particle 
by using a third-order transport matrix. In order to min- 
imize the multiple-scattering effects on the momentum 
resolution, the QQDQQ system was designed so as to 
make the < x\9 > term of the transport matrix to be 
zero between the focal planes close to the entrance and 
exit windows of 100 fim thick stainless steel of the vac- 
uum beam pipe and the drift chambers were made as 
thin as possible. The magnetic field of the dipole magnet 
(D2) was monitored throughout the experiment for every 
spill with a high-precision Hall probe in order to correct 
its fluctuation in the off-line analysis. 



±21 mm. The particle track with a position resolution 
of 300 ^m was measured by SDC's. 

In the present experiment, SKS was excited mainly 
to 272 A, which corresponds to 2.2 Tesla(T), in order 
to measure the spectrum in the quasi-free region cover- 
ing partly below the binding threshold. In addition, 
with Si and CII2 targets, SKS was excited to two other 
different settings of 320 A(2.4 T) and 210 A (1.9 T) to 
measure the spectrum well below the binding thresh- 
old and in the highly excited regions, respectively. The 
central momentum of SKS at 2.4 T, 2.2 T and 1.9 T 
were 780 MeV/c, 720 MeV/c and 620 MeV/c, respec- 
tively. The spectrum in one SKS setting could be con- 
nected smoothly to another as the acceptance regions 
partly overlap. The scattered-particle momentum was 
obtained particle-by-particle by reconstructing a particle 
trajectory with the Runge-Kutta tracking method using 
a precisely measured magnetic field map in each setting 
hj. The magnetic field was monitored with an NMR 
probe throughout the data acquisition in order to cor- 
rect its fluctuation in the off-line analysis. The fluctua- 
tion was as low as ±0.003%. 



B. SKS spectrometer 

The SKS spectrometer has a large acceptance of 100 
msr, a good momentum resolution of 0.1% (FWHM), and 
a linearity of within 0.1 MeV/c in its momentum accep- 
tance of ±10% The SKS system is thus suitable 
to obtain high statistics over a wide excitation energy 
region with a good energy resolution, while maintaining 
the sensitivity of the imaginary part of the S-nucleus po- 
tential because a bad resolution may cause the spectrum 
to smear out. 

SKS comprises a superconducting dipole magnet to- 
gether with four sets of drift chambers (SDCl-4) for mo- 
mentum reconstruction and three kinds of trigger coun- 
ters. The trigger counters comprise a scintillation counter 
wall (TOF), a Incite Cerenkov counter wall (LC) and two 
silica aerogel Cerenkov counters (ACl and AC2) as seen 
in Fig. n] The TOF counter comprised 15 vertical scin- 
tillation counters, each 7x100x3 cm'^ in size, and was 
used for scattered particle identification by measuring 
the time-of- flight from the reaction point. The LC was a 
threshold-type Cerenkov counter comprising 14 pieces of 
10x140x4 cm'^ Incite radiators (n=1.49), which discrim- 
inated protons from pions and kaons. ACl and AC2 were 
threshold-type silica aerogel Cerenkov counters (n=1.06) 
used to eliminate pions. The {'k~ ,K'^) trigger was de- 
fined as PIK= BEAM X TOF xLCx:4CT x AC2. The 
trigger rate was typically 300 counts for a beam rate of 
2.0x10^ per spill. 

The drift chambers SDCl and SDC2 were installed at 
the entrance of the SKS magnet, while SDC3 and SDC4 
were installed at the exit. SDCl and SDC2 had the same 
drift-cell structure as the BDC's, as they were exposed 
to the beam. SDC3 and SDC4 had a large drift space of 



TABLE I: Data summary of E438 



Target (s) 


SKS current 


Central momentum 


irradiated 




[Amp(TesIa)] 


of SKS [MeV/c] 


TT- [XIO''] 


CH2 and Si 


210 (1.9) 


630 


86.5 


CH2 and Si 


272 (2.2) 


720 


228.7 


CH2 and Si 


320 (2.4) 


780 


101.6 


CH2 and Ni 


272 (2.2) 


720 


252.6 


CH2 and In 


272 (2.2) 


720 


352.3 


CH2 and Bi 


272 (2.2) 


720 


225.2 


CH2 and Bi 


272 (2.2) 


720 


120.7 


CH2 only 


272 (2.2) 


720 


47.3 


CH2 only 


272 (2.2) 


720 


39.4 


Target empty 


272 (2.2) 


720 


94.0 



C. Experimental targets and Data summary 

In the present experiment we used CII2 (1.00±0.05 
g/cm^), Si (6.53±0.07 g/cm^), Ni (7.16±0.04 g/cm^). 
In (7.93±0.09 g/cm^) and Bi (9.74±0.10 g/cm^) of nat- 
ural isotopic composition as experimental targets. The 
quoted error on each target thickness came from the mea- 
surement. The thickness and the number of irradiated tt" 
on each target were optimized in order to obtain almost 
the same yield in the quasi-free region for production. 
The data were taken in two experimental cycles in 1999. 
Table n gives a data summary with specific SKS current 
settings according to the experimental requirement. The 
CH2 target was always put in tandem at a distance of 
250 mm upstream from other natural targets (Si or Ni 
or In or Bi), as demonstrated in Fig. El In order to cal- 
ibrate the horizontal axis (excitation energy scale) and 
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to check the rehabihty of the vertical axis (cross section) 
throughout the experiment, as well as to obtain the en- 
ergy resolution, the elementary, p{tt^ , K^)T,^ process, 
from the CH2 target was used. On the other hand, the C 
spectrum was also extracted from the CH2 target. The 
target empty data were used to check the background 
level in all inclusive spectra. 



IV. DATA ANALYSIS 

The hypernuclear mass {Mhy) in the {tt^ ,K^) reac- 
tion was calculated by the following relation: 



Mhy ^\I{E^+ Ma - Enf - {pI +Pk- 2p^PK cosO^k) 

(2) 

where E^^ and are the total energy and momentum of 
a pion; Ek and pk are those of a kaon. Ma is the mass 
of a target nucleus and 9 is the scattering angle of the 
reaction. The measured inclusive {'it~ ,K~^) spectra are 
presented as a function of the binding energy {B-^- ) of 
I]~ hyperon, which was defined as 

Ss- = Ma-1 + Ms- - Mhy , (3) 

where, Ma-i is the mass of a core nucleus at its ground 
state and M^- is the mass of a hyperon. 



A. (tt ,K'^) event selection and momentum 
reconstruction 

The main background in the PIK trigger was fast pro- 
tons that fired LC, while pions were well suppressed in 
the trigger level by ACl and AC2. At the first stage of 
the off-line analysis, a large background event in the PIK 
trigger was rejected by using only for the counter's infor- 
mation. An incident vr^ was selected from the time-of- 
fiight information between BHl and BH2, while the scat- 
tered was roughly selected from the ADC and TDC 
information of the TOP and LC counters. Then, the tt^ 
momentum and momentum were reconstructed from 
the track information of the BDC's and SDC's, respec- 
tively. In the tracking process, straight-line track candi- 
dates were first defined both at the entrance and exit of 
each spectrometer using a least-squares method. Then, 
the combination of the straight-line tracks that gave the 
least chi-square in the momentum reconstruction was as- 
signed as the best track candidate. In the beam spec- 
trometer, a third-order transport matrix was used for 
momentum reconstruction, while in the SKS spectrom- 
eter the momentum was calculated by the Runge-Kutta 
method with a measured magnetic field map. After mo- 
mentum reconstruction, the mass of a scattered particle 
was calculated as 
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FIG. 2: Schematic view of the target configuration. 
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FIG. 3: Typical distribution of the scattered-particles mass 
obtained from the CH2 and Si data. The gated region was 
used as good kaons. 
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FIG. 4: Z-vertex distribution after selecting good kaons from 
the CH2 and Si data. 



Mseat^^Vl^ , 



(4) 



5 



where /3 is the velocity of a scattered particle obtained 
from the time-of-flight and the flight path length between 
BH2 and TOF, and p is the reconstructed momentum. 
Figure 131 shows a typical mass spectrum of the scattered- 
particles obtained from the CH2 and Si data. Kaons were 
clearly identified separated from pions and protons, and 
the gated region was used as good events. 

The scattering angle and the vertex point were ob- 
tained from the local straight track in BDC3 and 4 for 
the incoming pion and a track obtained in the momen- 
tum reconstruction of the scattered particle in SKS. The 
Z-axis was defined as the beam direction. Figure0]shows 
a typical Z- vertex distribution from the CH2 and Si data 
obtained by selecting good kaons, where events scattered 
in both targets as well as in the BH2 counter were clearly 
identified. However, the Z-vertex resolution at very for- 
ward scattering angles was not good due to the multiple- 
scattering effects. A scattering angle [9) cut greater than 
4° was required to improve the Z-vertex resolution, as in 
Fig. El Two arrows indicated in Fig. 0] were the selected 
CH2 and Si events used in the analysis. 



B. Energy-scale calibration and energy resolution 

The calibration of the horizontal axis (binding energy, 
By,- ) was one of the very important points of the present 
experiment, because the quasi-free peak position of the 
inclusive spectrum is expected to provide information 
on the S-nucleus optical potential. Since the elemen- 
tary, p{Tr~ , K~^)'E~ process from the CH2 target was 
used for this purpose, the elementary peak stands at 
259.177±0.029 MeV/c^ which is the mass difference be- 
tween a S]~ and a proton |23|. The calibration process 
mainly consists of two parts (a detailed description can 
be found in Ref. 

First, the energy losses of the incident pions and out- 
going kaons in the target (s) were estimated from the 
beam through data taken with and without a target; sec- 
ond, correlations between the kaon momentum and the 
incident angles (horizontal and vertical) are needed to 
be solved so as to determine the SKS momentum offset 
value. This offset value should be kept the same for all 
sets of data in order to avoid any ambiguity in the cali- 
bration. 

Figure [SI shows a typical missing-mass spectrum from the 
CH2 of CH2 and Si data, where the kinematics was solved 
while considering the reaction with the proton target in 
CH2. The spectrum is plotted as a function of M^- - 
Mp in MeV/c^ at a scattering angle of Ok = 6° ± 2°, 
where a calibration of the horizontal axis was not done. 
As a result, the elementary peak position was not found 
at the expected position (259.177 MeV/c^). A small 
satellite peak (left) came from events which did not pass 
through the downstream Si target, as can be seen from 
Fig. [3 A broad bump in the spectrum corresponds to the 
events from C in the CH2 target. Two kinds of elemen- 
tary events were identified in order to take into account 



the energy losses correctly. Figure shows the horizon- 
tal axis corrected spectrum, where it was fitted with two 
Gaussians, as plotted with the solid curve. One Gaussian 
corresponds to the elementary peak, whereas the other 
one corresponds to the carbon contribution. From the fit- 
ting result, the elementary peak position was obtained at 
259.23 ± 0.13 MeV/c^, as expected, whereas the energy 
resolution was obtained to be 3.31 ± 0.33 MeV (FWHM). 
The elementary peak position and the energy resolution 
were checked for all data sets (Tabled, while keeping all 
of the offset parameters the same. A calibration of the 
horizontal axis was successfully achieved, as summarized 
in Table m where the elementary peak position was al- 
ways found to be the expected position within an error 
of ±0.2 MeV/c^. The error for the peak position in each 
data set came from the fitting result due to the statistical 
fluctuation. 

The linearity of the SKS momentum was checked by a 
7r+ beam through data without any target for several cen- 
tral beam momenta of 630 MeV/c to 830 MeV/c, where 
the SKS was excited at 272 A. Here, we assumed that the 
central momentum of the beam spectrometer was exactly 
proportional to the magnetic fields of the beam-line mag- 
nets. The linearity was found to be better than ±0.072 
MeV/c in the momentum range [2ll |. 



TABLE II: Elementary p{iT~ , K'^)'E~ peak position and the 
energy resolution after a horizontal axis calibration for all sets 
of data. The quoted errors are statistical. 



Data set and 
SKS current 



Elementary peak Resolution(FWHM) 
position(MeV/c^) (MeV) 



CH2 only (272A) 
CH2 and Si (210A) 
CH2 and Si (272A) 
CH2 and Ni (272A) 
CH2 and In (272A) 
CH2 and Bi (272A) 



259.20 ±0.10 
259.23 ±0.13 
259.27 ±0.14 
259.45 ±0.19 
259.39 ± 0.24 
259.62 ± 0.26 



2.06 ±0.20 

3.31 ±0.33 

3.32 ±0.29 
4.44 ±0.42 
4.79 ±0.50 
5.16 ±0.53 



C. Cross section 

In order to obtain the cross section from the experi- 
mental yield, the experimental efRciencies and SKS ac- 
ceptance were estimated. Then, the cross section was 
calculated as 



A 



K 



da 



dQ J ipx) ■ Na J2t It • 



eSDC12 ■ £SDC3i) • f decay ' dfl {p, 6) 



, (5) 



where A is the target mass number, px the thickness 
of the target in g/cm^, Avogadro's number, the 
yield number, li the number of irradiate tt", the sev- 
eral experimental and analysis efhciency factors: data- 
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FIG. 5: Missing-mass spectrum of the (vr" ,K^) reaction from 
CH2 in the proton-target kinematics before a horizontal-axis 
cahbration. See text for details. 
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FIG. 6; Horizontal-axis corrected missing-mass spectrum 
for the elementary process from CH2, where the peak po- 
sition is found at 259.23±0.13 MeV/c^ with a resolution of 
3.31±0.33 MeV (FWHM). The expected peak position is at 
259.177±0.029 MeV/c^ the mass difference between a E~ 
and a proton |20|l . 



acquisition efRciency, BDC's efficiency, K6 tracking ef- 
ficiency, beam-normalization factors, SKS tracking effi- 
ciency, Z-vertex cut efficiency, TOF and LC's efficiency, 
AG's accidental veto factor and TOF and LC's multiplic- 
ity cut efficiency. esDCi2 and esDCSi ai'e the efficiencies 
of SDC12 and SDC34, respectively. These chambers had 
an incident position dependence of the efRciency, which 
was estimated event by event, as explained below, /decay 
is the kaon decay rate and dn{p, 9) the effective solid 
angle of SKS as a function of the kaon momentum and 
scattering angle. 



1. Experimental and analysis efficiencies 

The efficiencies mentioned above were estimated sepa- 
rately for all data sets. However, typical efficiencies ob- 
tained from the CH2 and Si data with SKS at 272 A are 
described below. 

The beam- normalization factor (/beam) represents the 
fraction of 7r~ out of the total number of the beam 
{Nbeam), and was estimated from the BEAM trigger 
events as fbeam = (l-//i)-(l-/occ), where is the muon 
contamination in the 7r~ beam and face is the accidental 
coincidence rate between BHl and BH2. The e~ con- 
tamination was rejected by the good performance of cGC 
with an efficiency better than 99.9%. By requiring a good 
timing coincidence between BHl and BH2, other contam- 
ination in the pion beam was rejected. However, muons 
in the beam that were the decay products of pions, could 
not be separated from pions. We took this value from 
previous experiments, since it was studied and found to 
be 6.2% with a measurement error of ±2.0% for a 1.0 ^ 
1.10 7r+ beam [23, H^. Factor face was estimated as 



(Nb 



-N. 



BH1.2 



)/Ni 



beam : 



where Nbhi.2 is the number of 



events for which the time-of-flight between BHl and BH2 
was proper for 7r~. Factor face was typically 4.4±1.7%. 
Finally, the beam- normalization factor (/fceom) was found 
to be 89.6±2.6%. 

The BDC's efficiency was the total efficiency to obtain 
a straight track at the entrance and exit of the QQDQQ 
magnets, and was estimated to be 88.9±1.6%, whereas 
the K6 tracking efficiency was the analysis efficiency to 
reconstruct a particle trajectory in the beam spectrome- 
ter, which was found to be 96.5±1.4%. 

In the analysis, more than two hits on TOF or LC 
were rejected in order to reduce the background level. 
The corresponding efficiency for this cut was estimated 
to be 98.7±1.3%. The intrinsic efliciency of TOF and LC 
together was 99.6±0.3%. 

Due to the high counting rate of ACl and AC2, some 
of the PIK events were killed accidentally. The coinci- 
dence width between ACl, AC2 and BEAM x TOF xLC 
was 56±5 nsec, which was the dead time in the PIK trig- 
ger. The single counting rate of ACl and AC2 was typi- 
cally 1.5x10^ s~^. Thus, the AG's accidental veto factor 
{fAc) was calculated to be 99.1±0.2%. 

The SKS cut efficiency was estimated from the 
(7r~,p) events mixed in the PIK trigger. At first, a 
loose cut was applied to reconstruct particle trajecto- 
ries in SKS, and at the final stage of analysis a further 
tight cut was applied in order to reduce the background 
level as much as possible. The resolution in the parti- 
cle identification was significantly improved by this cut, 
where the corresponding cut efficiency was estimated to 
be 92.2±2.9%. The PID cut efficiency for kaons after the 
SKS cut was found to be 99.4±1.6%. 

The Z-vertex resolution at very forward scattering an- 
gles was bad due to the multiple-scattering effect, as also 
mentioned earlier. In addition to the scattering angle cut 
(>4°), tighter cuts were applied, as shown by the arrows 
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in Fig. 0] The corresponding cut efficiencies were ob- 
tained by fitting the spectrum with a Lorentz function. 
Those were 90.0±3.9% and 92.2±2.8% for CH2 and Si, 
respectively in the CH2 and Si data. 

The efficiency of SDCl and SDC2 [630012) was the 
total efficiency, including the analysis efficiency, and was 
estimated using the BEAM trigger. In SDCl and SDC2, 
the beam-counting rate per wire was quite high because 
the beam was focused at targets(s). As a result, a small 
degradation of the efficiency near the beam spot was ob- 
served. Therefore, the efficiency was estimated as a func- 
tion of the horizontal incident position for the event-by- 
event correction. The average efficiency of SDCl and 
SDC2 was typically 92.8±1.2%. 

The efficiency of SDC3 and SDC4, including the analy- 
sis efficiency, was estimated from the {tt~ ,p) events. Since 
it had an incident position dependence of the efficiency 
due to some noisy channels, the efficiency as a function of 
the horizontal incident position was used for corrections 
event by event, where the average value of the SDC3 and 
SDC4 efficiency was found to typically be 85.2±1.7%. 



TABLE III: Elementary p(7r", A'+)E" cross section from dif- 
ferent sets of data. The quoted errors are both statistical(stat) 
and systematic (sys). 



Data set and 
(SKS current) 



TT p — l-S ■ 

[fih/sr] 



CH2 only (272A) 
CH2 and Si (210A) 
CH2 and Si (272A) 
CH2 and Ni (272A) 
CH2 and In (272A) 
CH2 and Bi (272A) 



127.20 ± 9.59(stat) 
128.41 ± 8.13(stat) 
124.73 ± 7.48(stat) 
127.46 ± 6.52(stat) 
122.82 ± 5.42(stat) 
123.98 ± 5.14(stat 



±7.19(sys) 
±6.09(sys) 
±6.15(sys) 
±8.51(sys) 
±6.54(sys) 
±7.02(sys) 



from 6 — ^A9 to 6* -I- ^A6 in the polar angle, from to 
2(/) in the azimuthal angle, and from p — ^Ap to p+ ^Ap 
in the momentum. 



D. Elementary, p(7r ,K^)'E cross section 



2. Other factors 

The data-acquisition efficiency is caused by the dead 
time of the data-acquisition system. It was estimated 
from the ratio of the number of events accepted by the 
data-acquisition system to that of triggers, and was typ- 
ically 80.4±0.1%. 

The kaon decay rate in SKS was studied in detail with a 
Monte-Carlo simulation by T. Hasegawa l23|. The corre- 
sponding correction factor was typically 40.0±2.0%, and 
was corrected event by event while taking into account 
the momentum and the flight path length. 

3. SKS acceptance 

The acceptance of the SKS spectrometer to obtain the 
effective solid angle (dil) was calculated by a Monte- 
Carlo simulation code, GEANT The geometrical 
condition of the experiment, the effect of the energy loss, 
the multiple scattering effect and the off-line analysis cut 
condition were taken into account. In addition, a timing 
cut at the lower momentum side of the TOE wall was 
applied, as it was found by comparin g th e timing spec- 
tra from the simulation and the data |2l| . In the event 
generator, the distribution of the beam profile obtained 
from the experimental data was produced and the effec- 
tive solid angle was averaged on the distribution. It was 
calculated as a function of the scattering angle (6) and 
momentum (p), as follows: 

pe+(i/2)Ae I-2TV 
dn{p,d)^ dcos6 d,pxC{p,0), (6) 

Je~{i/2)Ae Jo 

where C (p, 9) is the ratio of the accepted events to the 
generated events. The events were generated uniformly 



Taking into account all of the above mentioned efficien- 
cies and the SKS acceptance, the differential cross section 
at the kaon scattering angle of 6° ± 2° in the laboratory 
frame was obtained as follows: 

V*^^ 740-8° Je=i° \dfl J J 8=4° 

Data from the elementary process in the present experi- 
ment was used mainly to calibrate the excitation energy 
scale as well as to check the accuracy of the measured 
cross section. The calibration of the energy scale has al- 
ready been described earlier. The elementary cross sec- 
tion obtained from all sets of data are presented in this 
section. In Fig. El the vertical axis is calibrated to be the 
cross section in units of /ib/sr/MeV/c^. By subtract- 
ing the carbon contribution, the elementary cross section 
was extracted from the proton peak and was found to 
be 128.41±8.13 ^b/sr/MeV/c^. Elementary cross sec- 
tions obtained from all sets of data are summarized in 
Table ITlll along with both the statistical and systematic 
errors. The elementary cross section was found to be con- 
sistent in each set of data within the errors. This demon- 
strates the stability of estimating the SKS acceptance as 
well as all of the detectors and analysis efficiencies during 
the experiment. 

The angular distribution of the elementary cross sec- 
tion was obtained over a wide angular range from the 
CH2 and Si data with SKS at 210A. Only this setting 
could well cover the kaon momentum over a wider an- 
gular range. In Fig. \7\ the obtained angular distribution 
is shown by empty circles, where the two cross points 
are that from old bubble-chamber data measured at a 
beam momentum of 1.225 GeV/c [2^. The past data 
was originally reported in the center-of-mass system. It 
was transfered to the Lab system in order to compare it 
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2 4 6 8 10 12 14 16 18 20 
Scattering angle (0|^+) (deg) 

FIG. 7: Angular distribution of the elementary reaction in the 
Lab system. The empty circles are the present measurement, 
where the crosses are from previous bubble-chamber data [2R| . 
See the text for details. The quoted errors are statistical. 



with the present result and two data points were found 
at 9K=iO° and 17.5°. The dotted curve was obtained 
by fitting the previous data with a Legendre function, 
where the sohd lines show the error boundaries. The 
present measurement agrees well with the previous one 
with better statistics. 



E. Systematic errors 

In order to obtain the total systematic errors of the 
cross section, the measurement error of the target thick- 
ness, the SKS acceptance, and the incident position de- 
pendent efficiencies of SDC12 and SDC34 were taken into 
account. The effects of other sources were found to be 
very small, so as to be negligible. The efficiencies of 
SDC12 and SDC34 varied with the position. Particularly, 
SDC34 had a significant inefficiency for some channels at 
the lower momentum side. In order to obtain systematic 
errors, the efficiency functions of SDC12 and SDC34 were 
varied within the fitting errors j21| . The SKS acceptance 
was also found to be affected at the lower momentum side 
due to the trigger timing problem, as also mentioned ear- 
lier. To study the systematic error, the SKS acceptance 
was recalculated by changing the TOF timing cut within 
the typical time resolution of the TOF counter. The to- 
tal systematic errors obtained for the elementary cross 
section are summarized in table Hill The systematic er- 
rors on inclusive spectra were found to be larger in the 
higher excitation energy region as compared to the lower 
excitation energy region, because the systematic errors 
were mainly affected at the lower momentum side. 



F. Background level 

One of the main advantages of the {n^ ,K^) reaction 
over the {K~,7t^) reactions is its low background na- 
ture. To examine the background level exactly, in ad- 
dition to the data on nuclear targets, some data were 
taken without any target, as shown in Table ^ where 
the SKS current setting was at 272 A in order to see the 
background level mostly near to the binding energy 
threshold. The target-empty data were analyzed using 
the same analysis program and conditions as that for the 
normal data. The background level was found to be very 
low, as can be seen in Fig.|Hl(broken line histogram). The 
solid line histogram was the {tt^ ,K'^) spectrum from the 
Si target. The horizontal axis is the binding energy of S~ 
{-B^-), whereas the vertical axis shows the counts. The 
total entries are 7655, and 26 for Si and target-empty 
data, respectively, where the number of tt" irradiations 
for target-empty data was about half the amount com- 
pared to that for the Si target. For a comparison, the 
vertical axis was normalized by the total number of the 
beam (iV^r- ) for each case. As can be seen, the back- 
ground was found to be almost uniform, about 2 orders 
of magnitude lower than the spectrum with the target, 
and there was no background event around the bound 
region. Thus, we neglected the background in analyzing 
the inclusive spectra for all of the targets. 



Entries 7655 




200 



-B;-- (MeV) 

FIG. 8: Background level in the {tt' ,K'^) reaction from the 
target-empty data analysis. The horizontal axis shows the 
binding energy of a E~, whereas the vertical axis shows the 
counts. The solid- line histogram is from the Si target, while 
broken-line is from the target-empty runs. 



V. EXPERIMENTAL RESULTS 

The inclusive {n^ ,K^) spectra at a scattering angle of 
6°±2° on Si, C, Ni, In and Bi are presented here. The 
differential cross section was obtained by Eq.[7| All of the 
inclusive spectra are also presented in the tabular form 
in Tables Ivml Hxl Ixl and IXll 
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A. Si spectrum 

Figure^lshows the inclusive {n^ ,K^) spectra on Si ob- 
tained from three data sets of CH2 and Si with three SKS 
current settings. The horizontal axis is the missing mass 
in terms of -B-^- , and the vertical axis is the differential 
cross section in /ib/sr/MeV. Three spectra from three set- 
tings were found to be matched very well for the overlap 
region of the acceptance. These three settings together 
cover a wide energy region so as to understand the gross 
feature of the whole spectral shape. By taking a weighted 
average from three spectra, the combined Si spectrum is 
shown in Fig. 1101 where the spectrum has been shown in 
both statistical and systematic errors. A small amount 
of CH2 contamination found in the Si spectrum was sub- 
tracted with a proper scale. The CH2 contamination was 
found in all inclusive spectra, and was subtracted by the 
same manner 0|. As can be seen in the figure, the cross 
section was found to gradually increase with an increase 
of -By,-, and the maximum was at around -By,- — 120 
MeV. There is also a significant yield found below the 
binding energy threshold. 



K* Momentum (MeV/c) 
850 790 730 670 605 535 




Binding energy (-Bj--) (MeV) 



FIG. 10; Si spectrum obtained as a weighted average from 
three spectra measured with three SKS current settings. The 
quoted errors are both statistical and systematic, where the 
boundaries of the statistical errors are shown by arms and the 
systematic errors are drawn beyond the statistical errors. 
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FIG. 9: Inclusive {-k' ,K^) spectrum on Si with three set- 
tings of the SKS current and with statistical errors only. The 
scattered kaon momentum is also depicted on the top of the 
figure. 
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FIG. 11: Inclusive (tt ,K ) spectrum on C from the CH2 
data. See the text for details. 



B. C spectrum 



the spectrum from a deeply bound region to around 80 
MeV above the binding threshold. 



The inclusive (7r",if+) spectrum on C was extracted 
from the CH2 data. Like the Si spectrum, the C spectrum 
was also taken with three SKS settings, and was found to 
be matched well. The averaged C spectrum is shown in 
Fig. ^3 The sudden gap in the spectrum at around 100 
MeV was due to removing a sharp peak in the CH2 data, 
which included both H and C events; C events were very 
hard to separate. Nevertheless, we can fairly well discuss 



C. Ni, In and Bi spectra 



Inclusive spectra on Ni, In, and Bi were taken only 
at an SKS current setting of 272 A, which are shown in 
Figs. El El and El respectively. 
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FIG. 12: Inclusive (tt spectrum on Ni. 
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FIG. 15: Comparison of all the inclusive [n^ ,K^) spectra on 
each target taken with an SKS current setting of 272A. The 
quoted errors are statistical one only. 
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FIG. 13: Inclusive (tt ,K'^) spectrum on In. 
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FIG. 14: Inclusive (tt ,K'^) spectrum on Bi. 



TABLE IV: Ratio of the cross section for each target as 
compared to the cross section on Si at an energy region of 
< —Be- < 90 MeV. The quoted errors are statistical. 



Data 



Ratio 



C 
Si 
Ni 
In 
Bi 



0.81±0.15 

1.0 
1.16±0.07 
1.26±0.09 
1.49±0.12 



D. Similarity of the spectra 



From an analysis of all the inclusive spectra, it was 
found that all spectra show a similarity in shape at least 
up to —By,- =90 MeV, as shown in Fig. ^| All spectra 
with only an SKS setting of 272A were used to make the 
comparison fair. 



1. Mass-number dependence of the cross section 



The mass- number dependence of the production 
cross section was obtained from the ratio of the inclu- 
sive Si spectrum to that of the others. The ratios were 
taken at a energy region of < —B^- < 90 MeV. The 
ratios were found to be very flat over the energy region, 
and were fitted by a straight line. The fitting results are 
summarized in Table HVl 
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potential for the S-nucleus potential, parameterized as 

fi2 1 d/ 



V{r) = (Fo + iWo)/(r) + yso 



dr 



f{r) = 



1 



(8) 

1 + cxp ((r ~ c)/z) ' ^ ^ 

where Vso and Vc denote the spin-orbit and Coulomb 
potentials, respectively. The wave function of the initial 
nuclear state was calculated with the Woods-Saxon type 
potential of the same parametrization. To avoid any con- 
fusion, "S" is superscripted to the parameters, like 
or W(f , if it is for the S-nucleus potential. 

A. (7r~,A'+) spectrum on Si 



FIG. 16: Mass-number dependence of the cross section com- 
pared with the eikonal approximation. The quoted errors are 
only statistical ones. 



2. Comparison of the mass-number dependence to the 
eikonal approximation 

The effective nucleon number [N^ff) of the ,K^) 
reaction on C, Si, Ni, In and Bi at 6° was calculated 
by the eikonal approximation. A brief description of the 
calculation is given in Ref. [2ll |. 

The mass-number dependence of the cross section ob- 
tained from the present data was compared to that calcu- 
lated by the eikonal approximation, as shown in Fig. 1161 
The horizontal axis is the target mass number, whereas 
the vertical axis is the ratio of the cross section as com- 
pared to the Si fTable UVjl . As can be seen in the figure, 
the mass-number dependence of the ,K^) reaction 
obtained from the present data shows a rather weak de- 
pendence compared to that found in a calculation by us- 
ing the eikonal approximation. The present data were 
fitted by a function of Ca x where Ca is a constant 
and a is the fitting parameter. The value of a was found 
to be 0.20±0.04. Because a can reflect the effect of dis- 
tortion, the present data show a stronger distortion than 
that from the eikonal approximation. 



VI. COMPARISON WITH THE DWIA 
CALCULATION 

We calculated the inclusive {■k~,K^) spectra within 
the framework of DWIA, and the calculated spectra were 
compared with the measured ones in order to obtain in- 
formation on the S-nucleus potential. The formalism of 
the calculation can be found in Appendix A. In the cal- 
culation, we assumed the Woods-Saxon type one-body 



We made a shape analysis of the {ir^ .K^) spectrum 
on Si, where the measured spectrum was fitted by the 
calculated ones for various S-nucleus potentials with the 
magnitude of each spectrum being a free parameter, as 
reported in Ref. [2^. We made an improvement on the 
calculation from those demonstrated in a previous re- 
port |26j| : we took into account a much higher angular 
momentum transfer, up to 21, in the reaction, by which 
the magnitude of the (tt", K^) spectrum at higher — -B^- 
was increased. This caused a change in the fitting result 
from that obtained in the previous report 26]. Figure lTTI 
shows the [tt^.K^) spectrum on Si fitted by the cal- 
culated spectra. All of the data points (57 points) that 
appear in the figure are employed for the fitting. The po- 
tential parameters used for the calculation are the same 
as those in (2^ , as listed in the second column of TablelVl 
Typical and dependences are shown in (a) and 
(b) in the figure. The spectrum shape shifts to a higher 
excited region with an increase of , and is weakly de- 
pendent on Wq . The slope of the calculated spectra with 
lower just above —B^-—0 is relatively steep, and 
thus does not reproduce the measured one. Particularly, 
the case of a shallow potential with (V^q^, W^) — {—10 
MeV,— 10 MeV), which were not excluded in the previous 
analysis in ^^C (stopped K~,7r^) |43|, is hard to repro- 
duce the measured spectrum. Figure [TBI shows a contour 
plot of the chi-square distribution. Here, the best-fit po- 
tential is found to be {V^ , W(f )=(110 MeV,-45 MeV), 
where the chi-square over the number of degree of free- 
dom is x^/DOF=29. 48/56, However, a wide region of Vip 
and Wq is still hardly excluded with respect to the confi- 
dence level. Therefore, the framework of the present the- 
oretical analysis only demonstrates that: (1) a strongly 
repulsive S-nucleus potential is required to reproduce the 
whole measured spectrum in shape, (2) while an attrac- 
tive potential is hard to reproduce the spectrum, and (3) 
the spectrum shape is not very sensitive to the imaginary 
potential, although a deep (~— 45 MeV) one seems to be 
favored. 

Since the strength of the wave function outside the nu- 
cleus plays an important role on the magnitude of the 
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FIG. 17: Inclusive 
calculated spectra. 



60 76 \QQ 

{tt~,K'^) spectrum on Si fitted by the 
Typical V(p and Wq" dependences are 
shown in (a) and (b). Curves in (a) are for Vo^=10~150 MeV 
every 20 MeV fixing at W(f =— 45 MeV except for the case 
of {V^, Wo )={-W MeV, -10 MeV). Curves in (b) are for 

Wo=-55 15 MeV every 10 MeV fixing at Fo^ = 110 MeV. 

The potential parameters used for the calculation are listed 
in the second column of Table [V] 




l^f (MeV) 



FIG. 18: Chi-square distribution in Vq^ and Wq for Si. The 
potential parameters used for the calculation are listed in the 
second column of Table [V] 



TABLE V: Potential parameters used for the DWIA calcula- 
tion on the (tt" , A""'')reaction. 



Parameter\target 


2«Si 




"^In 


209-Q- 


Us" 










Vo{MeY) 


-10~4-150 




+90 




Wo (MeV) 


-60~-10 




-40 




Vso(MeV) 
c(fm)'' 











3.3 


4.23 


5.33 


6.52 


2(fm) 




0.67 






Ut" 










Vo(MeV)'' 


-49.6 -54.4 


-51.6 


-51.4 


-55.5 


Wo (MeV) 











Vso(MeV) 




7 






c(fm)'' 


4.09 3.82 


4.95 


6.24 


7.42 




0.536 


0.517 


0.563 


0.468 



"S-nucleus potential. 
''c=l.lx(A- 1)1/3. 

''proton single-particle potential of the initial nucleus. 

''these are adjusted to repr odu ce the proton separation energy l27l 
and nuclear radius (r^)l/^ l28l . For Si, a smaller potential radius 
taken to be 3.82 fm was also calculated for comparison. 
■^Taken from 



inclusive spectrum in the repulsive S-nucleus potential, 
the potential radius is expected to be sensitive to the 
spectrum shape. Thus, we made another set of calcula- 
tions with a radius parameter of c=3.82 fm and a depth 
of Vb=— 54.5 MeV in the initial Si potential, as listed 
in the third column of Table Ivl The parameter c=3.82 
fm was chosen to be c«1.26xyl^/^, which is close to the 
usually quoted radius parameter of the optical model. 
Fitted spectra are shown in typical cases of and Wq 
in Fig. El Figure EOl shows a contour plot of the chi- 
square distribution. In this case, the best-fit potential 
is found to be {V^, Wo^)=(90 MeV, -40 MeV), where 



X^/DOF=27. 33/56. In Fig. [201 a good chi-square region 
becomes lower in the axis, compared with Fig. E| 
With this choice of the potential parameters, we again 
demonstrate that a repulsive potential is required to re- 
produce the measured spectrum in shape, while an at- 
tractive potential is still hard to reproduce the spectrum. 
The spectrum shape is weakly sensitive to W^. However, 
this result suggests that the choice of the potential radius 
parameter causes a systematic shift of the favored region 
in V;^ and . 
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FIG. 19: Same as Fig. 1171 but the curves in (a) are for 
Vo^ = 10~150 MeV every 20 MeV fixing at Wo=-40 MeV, 
except for the case of {V(p, W'o^)=(-10 MeV, -10 MeV). 
Curves in (b) are for iyo^ = -60~-10 MeV every 10 MeV fix- 
ing at V(p=90 MeV. The potential parameters used for the 
calculation are listed in the third column of Table [V] 



B. Ni, In, and Bi 

The fitting result in Si was applied to the other heav- 
ier targets. Figure ET1 shows the calculated spectra with 
(V'o^, Wo^)=(90 MeV,-40 MeV) fitted to the measured 
ones in Ni, In, and Bi. The other potential parameters 
are listed in Table Ivl The magnitude of the spectrum is 
arbitrarily adjusted, and all of the plotted data points in 
the figures were employed for the fitting. The values of 
xVDOF in Ni, In, and Bi are 54.65/52, 33.51/53, and 
37.55/48, respectively. The results show that the calcu- 
lated spectra do not contradict the measured ones. We 
also applied a chi-square test, as was done for Si. The fit- 




(MeV) 



FIG. 20: Same as Fig. 1181 but the potential parameters used 
for the calculation are listed in the third column of Table IVl 



ting results show rather more repulsive potentials at the 
best chi-square. However, it is unclear if the measured 
spectra maintain a sufficient sensitivity for a very repul- 
sive potential greater than lg^~130 MeV, or so, since 
the data quality may decrease above -Bi;~130 MeV due 
to the limited acceptance of the spectrometer, as dis- 
cussed in Section V. It is not very meaningful to dis- 
cuss the importance of the better chi-squares compared 
to those obtained in cases of {V^, W^)={90 MeV,-40 
MeV) because the test only check the inconsistency 
of the model. Therefore, we only demonstrate here that 
the repulsive potential obtained from the fitting analysis 
in Si is also applicable to the cases of Ni, In, and Bi. 



C. (vr^jif^) spectrum on C 

We also applied the above-mentioned formalism to the 
inclusive (tt^, K'^) spectrum on ^^C at the same incident 
pion momentum of 1.2 GeV/c. In this case, the spectrum 
is dominated by the quasi-free A production. Since the 
A-nucleus potential is well known, it is a good test to 
check if this formalism is applicable or not. The (tt"*", iir+) 
spectrum on ^^C is plotted together with the calculated 
spectra in Fig. The potential parameters used for the 
calculation are summarized in Table lVll The depth of the 
A-nucleus potential was taken to be yo'^=-30 MeV (solid 
line). The spectrum for V^'^=0 MeV was also calculated 
for a comparison (dashed line). Many lambda hypernu- 
clear states exist at around the lambda binding threshold 
and form a complicated structure in the spectral shape. 
More careful calculations on the relevant wave functions 
of the states, the width of the states, Gaussian convolu- 
tion to simulate the experimental resolution (Lorcntzian 
convolution in the present framework), and so on would 
be required to reproduce the spectrum much better in 
this region. Fractions of the S'''"'" productions contribute 
to the spectrum above — i?A~80 MeV. For this reason, we 
avoided to apply the x"-fitting test. The measured spec- 
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FIG. 21: Inclusive (vr , A'"'') spectra on Ni, In, and Bi fitted by 
the calculated spectra with {Vi^ ,Wo)^{90 MeV, -40 MeV). 
The other potential parameters used for the calculations are 
listed in the fourth to the sixth columns of Table IVl 



trum is thus overlaid with the calculated ones, where the 
magnitude of which are arbitrarily adjusted by eye. The 
solid line reproduces the gross feature of the spectrum 
over a wide range of the A-binding energy. 



VII. DISCUSSION 

The measured spectra are characterized by gradually 
increasing the cross section along with an increase of 
—By,- and the yield maximum at around —B^- as high 
as 120 MeV, as mentioned in Section V. Within the 
present framework of the calculation, we obtain a re- 
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FIG. 22: Inclusive (tt"*", A"*") spectrum on C is shown. The 
solid and dashed lines are calculated spectra with Vo^=— 30 
MeV and MeV, respectively. 



TABLE VI: Potential parameters used for the DWIA calcu- 



lation on the 


(tt"*", ii'+)reaction. 








A-nucleus pot. 


initial C 




Ua 




Ut'' 


VoiMeY) 


-30 





-50.3' 


Wo (MeV) 




Qd 


0" 


Vso(MeV) 


1 


1 


7 


c(fm) 


2.45^' 


2.45-'' 


3.15" 


2(fm) 


0.6 


0.6 


0.51 



"U^ was calculated for a comparison. 

'neutron single-particle potential of the initial nucleus (C). 

■^these are adjusted to reproduce the neutron separati on energy 
(18.731 MeV) |0 and nuclear radius (r2)l/2=2.454 fm I2I . 

''A constant value of —1.5 MeV was added to the imaginary po- 
tential to simulate the experimental energy resolution 

•^A constant value of —6 MeV was added to the imaginary poten- 
tial for the s-nucleon hole state. 

•'c=l.lx(yl- 1)1/3. 



pulsive E-nucleus potential to reproduce the measured 
spectra. A repulsive potential reduces the cross section 
at around the E^-binding threshold, and shifts the yield 
to a higher —B^- region. In the present analysis, we 
took into account only the leading term in the inclu- 
sive {tt~,K^) reaction. One can consider a multiple- 
scattering (multi-step) process in the 7r~ incident, 
production reaction, which is expected to increase the 
yield at a higher —B-^- . Although the effect of this multi- 
step process is open for future studies, it might reduce the 
size of the repulsive potential to fit the measured spectra. 

If the E-nucleus potential is repulsive, it is natural to 
interpret that yields in the E-bound region are due to 
a non-zero size of the imaginary potential, since these 
yields are beyond those expected from the Coulomb po- 
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tential and those mixed in with by smearing with the ex- 
perimental resolution. It is suggested that the imaginary 
potential is estimated to be Wq^~12 MeV or more by 
the standard expression Wq^T /2^—vaavPp{0)/2, where 
(vcav) denotes the Fermi averaging of the product of the 
relative velocity in Y.N and the T^N AN conversion 
cross section, and Pp(0) is the proton density at the nu- 
clear center [s^l- This estimation is rather good for a 
low-energy EiV interaction. For a higher momentum 
Y,N interaction, as is expected in a large momentum- 
transfer reaction, (tt" , K~^) , the imaginary potential may 
be deeper according to theoretical estimations |^ l32| . 
although this statement is still unclear, since the experi- 
mental data on the TjN AN cross section are limited 
at a higher momentum region. 

Several sets of the hyperon-nucleon potential, based on 
the one-boson-exchange model, have been reported from 
the Nijmegen group j29j. There have been several at- 
tempts to calculate the E potential in nuclei and/or in nu- 
clear matter using this two-body potential. Most of them 
gave rather attractive S potentials '3l','33','34','3^,'3^, but 
those with the Nijmegen model-F (NF) potential gave 
repulsive E-nucleus potentials [s^. l35l | . The size of the 
repulsive potential was estimated to be -1-3. 6-8. li MeV 
in nuclear matter by using NF 35] . They have reported 
the S-nucleus folding potentials in ^®Si-S, based on the 
Nijmegen model-D (ND), F, and Softcore (NSC) YN po- 
tentials • Among them, the ND-based potential is 
the most attractive, and the NSC-based one is weakly 
attractive. The NF-based potential has a repulsive core 
with a height ~30 MeV and an attractive pocket at the 
nuclear surface. The sizes of the imaginary potentials 

for ND, NF, and NSC are 3.5 MeV, -7.4 MeV, 

and —15 MeV. The E-nucleus potential suggested from 
the present measurement and the fitting analysis is more 
repulsive than those based on ND, NF, and NSC. A E 
single-particle potential has also been calculated us- 
ing the quark-model-based hyperon-nucleon interaction, 
model-FSS proposed by the Kyoto-Niigata group [37l |. 
The E potential turns out to be 20.4 MeV repulsive, since 
the SA^(/=3/2) '^Si state shows a strongly repulsive na- 
ture due to the Pauli-exclusion effect between the quarks. 

A repulsive E-nucleus potential has been derived so 
as to reproduce the S]~-atomic X-ray data by employ- 
ing ajjhenomenological density dependent (DD) poten- 
tial |38, 39], or the relativistic mean field (RMF) theory 
[3^. The DD potential shows a strongly repulsive core 
having a height of ~95 MeV at the nuclear center with a 
shallow attractive tail outside of the nucleus; the imag- 
inary potential is as deep as W(f ~-35 MeV. The RMF 
approach was used to construct a Schrodinger-equivalent 
E-nucleus potential from the scalar (attractive) and vec- 
tor (repulsive) Dirac potentials, while fitting the experi- 
mental E~-atomic X-ray shift and width. This approach 
shows a change of the E-nucleus potential from attractive 
to repulsive in the nuclear interior along with an increase 
of a^^ from 1/3 to 1, where a^^ is the vector meson cou- 
pling ratio. The fitting is better when a^j increases. In 
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FIG. 23: Inclusive {■k~,K'^) spectrum on Si fitted by the 
calculated spectrum with the DD potetnial. 



the case of Si, the real and imaginary potentials are, re- 
spectively, about 30~40 MeV, and as deep as 40 MeV in 
the nuclear interior. The E-nucleus potential derived by 
fitting the E~-atomic data seems qualitatively close to 
that obtained from the present analysis of the (tt", if+) 
spectra. However, the calculated Si{TT^,K'^) spectrum 
with the DD potential within the present framework did 
not fit the measured spectrum very well, as shown in 
Fig. 123 This does not show that the volume of the DD 
potential is strongly repulsive. Since the X-ray data sug- 
gests the attactive potential in the atomic orbital region, 
the present Woods-Saxon type repulsive potential does 
not reproduce the X-ray data. Therefore, further studies 
on the E-nucleus potential would be required to explain 
the present measurements and the X-ray data, where 
more careful choices of the nuclear and E-nuclear radii 
and the figure of the E-nucleus potential at the nuclear 
surface and outside the nucleus would be necessary. 

The hyperon constituents of neutron star cores are 
closely related to the hyperon single-particle potential 
in dense nuclear matter, and hence are discussed based 
on the results from hypernuclear studies 0, 0, |3^ ^3 ■ 
However, discussions concerning the hyperons other than 
A are ambiguous, particularly for E, because the interac- 
tions of the hyperon to nucleon and other hyperons are 
not well established. The E~ appearance in dense neu- 
tron matter is sometimes discussed based on an attractive 
E-nucle us p otential in nuclear matter of normal density 
[T3L lisl |40| ]. On the other hand, the hyperon mixing in 
dense neutron matter is discussed based on the repul- 
sive E-nucle us p otential derived from the E~ -atomic X- 
ray data in [iSj. In this regard, the repulsive E-nucleus 
potential from the present result should be taken into 
account for discussing the hyperon dynamics in neutron- 
star cores. 
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VIII. SUMMARY AND CONCLUSIONS 

In the present experiment, we measured the inclusive 
{t:~ ,K'^) spectra for the first time on several medium- 
to-heavy nuclear targets (C, Si, Ni, In, and Bi) with 
reasonable statistics. The calibration of the horizontal 
axis (-Be-) was done successfully with a precision of bet- 
ter than ±0.1 MeV; also the angular distribution of the 
measured elementary cross section well reproduced previ- 
ous bubble-chamber data, indicating a reliability of the 
measured cross section. The energy resolution was ob- 
tained to be 3.31±0.33 MeV (FWHM). All spectra show 
a similarity in shape, at least at an energy region of 
< —Be- < 90 MeV; also the mass-number dependence 
of the cross section was found to be rather weaker than 
that from the eikonal approximation. 

We perfirmed a DWIA calculation of the inclusive 
(tt" , spectra on Si, Ni. In, and Bi. for various depths 
(heights) of the Woods-Saxon-type, one-body S-nucleus 
potential, and compared it with the measured ones con- 
cerning the shape. The present analysis demonstrated 
that a repulsive S-nucleus potential having a nonzero 
size of the imaginary part is required to reproduce the 
observed spectra. However, further studies on the E- 
nucleus potential, particularly choices of the nuclear and 
S-miclcar radii and the details of the potential shape at 
the nuclear surface and outside the nucleus, would be re- 
quired in order to explain both the present measurements 
and the X-ray data. 
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TABLE VII: Inclusive (tt ,K+) spectra on Si as a table. The 
quoted values of -B^- are at the center of bins. 



-Bj^- o"4o_so Statistical Systematic 
(MeV) (/ib/sr/McV) errors errors 



-50.0 





.0308 





.0127 


0, 
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-46.0 


0, 
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0, 


.0109 
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0, 
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0, 
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0, 


.0576 


0, 
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0. 
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0, 
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0, 
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0. 
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0, 
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0, 
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0. 
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0, 
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0, 
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0. 
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0, 
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0, 
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0. 
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0, 
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0. 
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0, 
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0, 
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0. 
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0, 
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0, 
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0. 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 
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0, 


.0140 
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1 
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0, 
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1 


.0527 
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0, 
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1 
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0, 
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1 
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.0594 


0, 
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1 
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.0719 


0, 
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1 
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.0629 


0, 
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0, 
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0, 
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0, 
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0, 
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.0842 
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0, 
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0, 
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0, 
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0, 
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\ 
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0, 
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0, 


.1171 


0. 
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0, 
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0. 
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0. 
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0. 
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1, 
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174.0 


1, 
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0, 


.2263 
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0, 


.2822 


206.0 


0, 
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0, 


.1489 


0, 


.7257 


210.0 


0, 


.3247 


0, 


.2302 


0. 


.1518 
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TABLE VIII: Inclusive [ir' ,K+) spectra on C as a table. The 

quoted values of —B^- arc at the center of bins. Data in 
the region 82 MeV<— B^- <106 MeV are omitted because of 
contamination by the events from H of the CH2 target. 



-Bj^- cr4o_go Statistical Systematic 
(MeV) (;nb/sr/MeV) errors errors 
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TABLE IX: Inclusive (tt ,K+) spectra on Ni as a table. The 
quoted values of -B^- are at the center of bins. 



-Bj^- CT40-S0 Statistical Systematic 
(MeV) (/ib/sr/MeV) errors errors 
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0. 


.0009 


-6.00 





.0824 





.0234 


0, 


.0010 


-2.00 


0, 


.0924 


0, 


.0226 


0, 


.0012 


2.000 


0, 


.0837 


0, 


.0225 


0, 


.0009 


6.000 


0, 


.0722 


0, 


.0205 


0, 


.0009 


10.00 


0, 


.1197 





.0238 


0, 


.0014 


14.00 


0, 


.1442 


0, 


.0248 


0, 


.0016 


18.00 


0, 


.1342 





.0240 


0, 


.0016 


22.00 





.1588 





.0248 


0, 


.0016 


26.00 





.2233 


0, 


.0283 


0, 


.0025 


30.00 





.1985 


0, 


.0268 


0, 


.0021 


34.00 





.2382 





.0292 


0, 


.0026 


38.00 





.3454 





.0333 


0, 


.0041 


42.00 





.3787 


0, 


.0347 


0, 


.0049 


46.00 





.4000 


0, 


.0355 


0, 


.0047 


50.00 


0, 


.4529 


0, 


.0383 


0, 


.0077 


54.00 


0, 


.5349 


0, 


.0407 


0, 


.0099 


58.00 


0, 


.6191 


0, 


.0441 


0, 


.0154 


62.00 


0, 


.6318 


0, 


.0441 


0, 


.0181 


66.00 





.7529 





.0496 


0, 


.0290 


70.00 





.8350 





.0500 


0, 


.0355 


74.00 


0, 


.8329 


0, 


.0522 


0. 


.0403 


78.00 


0, 


.9241 


0, 


.0609 


0. 


.0556 


82.00 


1, 


.1104 


0, 


.0599 


0. 


.0686 


86.00 


1, 


.2185 


0, 


.0635 


0. 


.0821 


90.00 


1, 


.2426 


0, 


.0753 


0. 


.0761 


94.00 


1, 


.3347 


0, 


.0689 


0. 


.0654 


98.00 


1, 


.4305 


0, 


.0841 


0. 


.1260 


102.0 


1, 


.4465 


0, 


.0887 


0. 


.1298 


106.0 


1, 


.3691 


0, 


.0923 


0. 


.1248 


110.0 


1, 


.4438 





.0879 


0, 


.1044 


114.0 


1, 


.4801 


0, 


.1038 


0, 


.1194 


118.0 


1, 


.4224 


0, 


.0983 


0, 


.1876 


122.0 


1 


.3406 





.1075 


0, 


.1183 


126.0 


1, 


.4306 





.1707 


0, 


.1561 


130.0 


1, 


.3283 





.1895 


0, 


.1368 


134.0 


1, 


.5228 





.1890 


0, 


.2007 


138.0 


1, 


.3421 


0, 


.1820 


0, 


.2400 


142.0 


1, 


.2459 


0, 


.1937 


0, 


.0881 


146.0 


1, 


.1056 


0, 


.2547 


0, 


.2580 


150.0 


1 


.1761 


0, 


.3156 


0, 


.3438 


154.0 


1 


.2218 


0, 


.5745 


0, 


.4022 


158.0 





.6356 


0, 


.3673 


0, 


.2316 


162.0 





.9136 


0, 


.4570 


0, 


.1678 


166.0 


0, 


.4536 





.4538 


0, 


.2631 


170.0 


0, 


.9689 


0, 


.7014 


0. 


.2788 
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TABLE X: Inclusive (tt ,K'^) spectra on In as a table. The 
quoted values of -B^- are at the center of bins. 



-Bj^- CT40-80 Statistical Systematic 
(McV) (/ib/sr/McV) errors errors 



-50.0 





.0855 





.0353 


0, 


.0039 


-46.0 





.0587 


0, 


.0281 


0, 


.0017 


-42.0 


0, 


.0035 


0, 


.0143 


0. 


.0002 


-38.0 


0, 


.0499 


0, 


.0248 


0. 


.0009 


-34.0 


0, 


.0302 


0, 


.0226 


0. 


.0022 


-30.0 


0, 


.0415 


0, 


.0202 


0. 


.0021 


-26.0 


0, 


.0329 


0, 


.0194 


0. 


.0010 


-22.0 


0, 


.0490 


0, 


.0216 


0. 


.0021 


-18.0 


0, 


.0814 


0, 


.0258 


0. 


.0121 


-14.0 


0, 


.0724 


0, 


.0233 


0. 


.0030 


-10.0 


0, 


.0763 


0, 


.0233 


0. 


.0058 


-6.00 


0, 


.0860 





.0234 


0, 


.0035 


-2.00 


0, 


.0650 





.0224 


0, 


.0014 


2.000 





.0865 


0, 


.0247 


0, 


.0109 


6.000 





.0846 





.0224 


0, 


.0055 


10.00 





.1312 


0, 


.0275 


0, 


.0176 


14.00 


0, 


.1349 


0, 


.0235 


0, 


.0097 


18.00 





.1546 





.0264 


0, 


.0168 


22.00 





.1861 





.0258 


0, 


.0062 


26.00 





.2204 


0, 


.0327 


0, 


.0097 


30.00 





.2731 


0, 


.0356 


0, 


.0168 


34.00 





.3003 


0, 


.0370 


0, 


.0062 


38.00 





.3572 


0, 


.0412 


0, 


.0063 


42.00 


0, 


.4032 


0, 


.0442 


0, 


.0213 


46.00 


0, 


.3974 


0, 


.0416 


0, 


.0091 


50.00 


0, 


.5141 


0, 


.0476 


0, 


.0212 


54.00 


0, 


.6095 


0, 


.0521 


0, 


.0124 


58.00 





.7018 


0, 


.0558 


0, 


.0187 


62.00 





.6084 





.0516 


0, 


.0155 


66.00 





.7704 





.0586 


0, 


,0261 


70.00 





.8492 





.0645 


0, 


.0564 


74.00 


0, 


.8939 





.0662 


0. 


.0524 


78.00 


0, 


.9444 


0, 


.0669 


0. 


.0824 


82.00 


1, 


.1461 


0, 


.0755 


0. 


.0635 


86.00 


1, 


.2693 





.0829 


0. 


.1165 


90.00 


1, 


.3236 


0, 


.0901 


0. 


.1501 


94.00 


1, 


.3309 


0, 


.0924 


0. 


.1541 


98.00 


1, 


.5872 


0, 


.1001 


0. 


.1706 


102.0 


1, 


.6932 


0, 


.1170 


0. 


.1890 


106.0 


1, 


.6785 


0, 


.1125 


0. 


.1613 


110.0 


1 


.7317 





.1204 


0, 


.1089 


114.0 


1 


.4510 


0, 


.1092 


0, 


.1941 


118.0 


1 


.8136 


0, 


.1322 


0, 


.1464 


122.0 


1 


.8225 





.1271 


0, 


.1276 


126.0 


1, 


.7817 





.1504 


0, 


.1325 


130.0 


1, 


.7700 





.1312 


0, 


.1613 


134.0 


1, 


.4876 





.1593 


0, 


.3185 


138.0 


1, 


.6231 





.1696 


0, 


.2132 


142.0 


1, 


.8433 





.2599 


0, 


.2551 


146.0 


1 


.4337 





.2636 


0, 


,1433 


150.0 


1, 


.1025 





.2401 


0, 


,4085 


154.0 


2, 


.1866 


0, 


.5857 


0, 


,3323 


158.0 


1, 


.4880 





.5633 


0, 


,2154 


162.0 


1, 


.4392 





.7206 


0, 


,5350 


166.0 


3 


.8892 


2, 


.8229 


0, 


,9892 


170.0 


0, 


.9020 


0, 


.2712 


0. 


,2883 



TABLE XI: Inclusive {n' ,K+) spectra on Bi as a table. The 
quoted values of -B^- are at the center of bins. 



-Bj^- CT40-S0 Statistical Systematic 
(MeV) (/ib/sr/MeV) errors errors 



-50.0 





,0317 





,0257 


0, 


,0318 


-46.0 


0, 


,0692 


0, 


,0383 


0, 


,0120 


-42.0 


0, 


.0440 


0, 


.0260 


0, 


,0010 


-38.0 


0, 


.0262 


0, 


.0214 


0, 


,0005 


-34.0 


0, 


.0409 


0, 


.0300 


0, 


,0002 


-30.0 


0, 


.0353 


0, 


.0195 


0, 


,0009 


-26.0 


0, 


.0817 


0, 


.0339 


0, 


,0017 


-22.0 


0, 


.0511 


0, 


.0238 


0, 


,0129 


-18.0 


0, 


.0733 


0, 


.0295 


0, 


,0011 


-14.0 


0, 


.0787 


0, 


.0256 


0, 


,0021 


-10.0 


0, 


.0696 


0, 


.0236 


0, 


,0291 


-6.00 





.0899 





.0317 


0, 


,0154 


-2.00 


0, 


.0698 


0, 


.0247 


0, 


,0032 


2.000 


0, 


.0680 


0, 


.0243 


0, 


,0086 


6.000 


0, 


.1402 





.0296 


0, 


,0030 


10.00 


0, 


.1323 





.0283 


0, 


,0029 


14.00 


0, 


.1472 


0, 


.0292 


0, 


,0108 


18.00 





.1893 





.0325 


0, 


,0033 


22.00 





.2246 





.0385 


0, 


,0187 


26.00 





.2397 





.0381 


0, 


,0112 


30.00 





.2858 





.0397 


0, 


,0050 


34.00 





.3895 





.0473 


0, 


,0114 


38.00 





.4681 





.0538 


0, 


,0214 


42.00 





.4490 


0, 


.0521 


0, 


,0250 


46.00 





.5479 


0, 


.0567 


0, 


,0312 


50.00 


0, 


.5748 


0, 


.0595 


0, 


,0378 


54.00 


0, 


.6558 


0, 


.0635 


0, 


,0258 


58.00 


0, 


.6803 


0, 


.0660 


0, 


,0222 


62.00 


0, 


.9564 


0, 


.0791 


0, 


,0756 


66.00 





.9413 





.0802 


0, 


,1025 


70.00 


1 


.0287 





.0865 


0, 


,0360 


74.00 


1, 


.0385 


0, 


.0868 


0, 


,0627 


78.00 


1, 


.3791 


0, 


.1059 


0, 


,1103 


82.00 


1, 


.3769 


0, 


.1088 


0, 


,1643 


86.00 


1, 


.4988 


0, 


.1157 


0, 


,1420 


90.00 


1, 


.7458 


0, 


.1318 


0, 


,2109 


94.00 


1, 


.6178 


0, 


.1326 


0, 


,2023 


98.00 


2, 


.1074 


0, 


.1728 


0, 


,2810 


102.0 


2, 


.1719 


0, 


.1984 


0, 


,1371 


106.0 


2, 


.1885 


0, 


.2283 


0, 


,3826 


110.0 


1, 


.6772 





.1956 


0, 


,4280 


114.0 


2, 


.1402 


0, 


.1992 


0, 


,2697 


118.0 


2, 


.2187 


0, 


.2057 


0, 


,0944 


122.0 


2, 


.1832 





.2111 


0, 


,1191 


126.0 


2, 


.3295 


0, 


.2213 


0, 


,2073 


130.0 


2, 


.1486 


0, 


.2420 


0, 


,4238 


134.0 


2, 


.7211 


0, 


.2151 


0, 


,5652 


138.0 


1 


.9433 





.3171 


0, 


,4770 


142.0 


2 


.0725 





.3910 


0, 


,3771 


146.0 


2, 


.0192 





.6665 


0, 


,4241 


150.0 


2, 


.3897 


0, 


.8487 


0, 


,2567 


154.0 


1, 


.8168 


0, 


.7440 


1, 


,4051 



APPENDIX A: DWIA CALCULATION 
FORMALISM 



A distorted-wave impulse approximation has often 

been used to evaluate the cross sections of the hyper- 
nuclear states in {K~,tt^), {Tr'^,K^), and others. The 
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Morimatsu-Yazaki formalism has been successfully ap- 
plied to an analysis of the {K~,Tr^) spectra in helium 
|4lj | and carbon ji^. It has also been applied to the 
(tt", K'^) reaction. The cross section of the inclusive re- 
action can be written as 



(Pa 

dQ.dE 



da 
dVlpi 



SiE), 



(Al) 



where /3 is the kincmatical factor for a coordinate transfer 
from a two-body system to a many-body system |45| , and 
da/dQei represents the averaged differential cross section 
of the elementary reaction. Here, the elementary cross 
section {da/dfl{s,nK)) is averaged over the momentum 
of a proton (fc) moving in the target nucleus with a weight 
of the momentum distribution (p(fc)), as written as 



da 
dQfi 



(E) 



jpik)^{s,nK)s{k-p)dk 



J p{k)S{k ~ P)dk 



(A2) 



P = k 



K+ 



The delta function is required for energy-momentum con- 
servation in the elementary reaction. Relevant particles 
are treated on mass-shell, and the possible k are given 
by fc^- and kji+ at E{—B^~). As a result, da/dflei 
depends on E. Figure shows calculated da/dVlei as 
a function of —B^,- in Si. In this calculation, the ele- 
mentary [t:^ ,K~^) reaction cross section, da/dQ{s,D,K), 
was taken from Refs. H^H^], where s is the Mandclstam 
variable. The momentum distribution of a bound nu- 
cleon p{k) was given by the single-particle wave function 
of the hole state. For the elementary {tt~^,K^) reaction, 
da/dfl{s,flK) was calculated by using the parameters 
obtained from the partial wave analysis The ob- 

tained da/dflei in Fig. |31 reflects that the elementary 
(7r~,_ftr+) cross section increases when s is smaller at a 
forward scattering angle in the Lab. frame. This en- 
ergy dependence is important to explain the spectrum 
shape. Particularly, the (t:^ , K~^) spectrum shape can- 
not be reproduced without taking this energy dependence 
into account. 

The strength function {S{E)) is written as 

SiE)^ 

/ drdr'fUr)Gc.^,{E;r',r)U{r'), (A4) 



Ur) = x^->iR)x'^+HR)H^Nir)\i), 
R = {MjMHy)r. 



(A5) 
(A6) 



Here, fa is the form factor characterized by the distorted 
waves of the incident and outgoing particles {x^^^ a-nd 
X^~^) and the nucleon-hole state ((a|V'Ar(r)|i)). The re- 
coil effect of the residual nucleus is taken into consid- 
eration in Eq. IA6I The response function can be de- 
scribed by means of the Green's function (c./. eq. (2.2) 
in Ref. |4y|). Here, we assume the Woods-Saxon type 
one-body S-nucleus potential {U s{r)=V^ {r)+iW^ {r)) . 



400 



350 




-25 25 50 75 100 125 150 175 
FIG. 24: Fermi- averaged elementary p{'K~ , K^)Ti~ cross sec- 



(A3) 

tion, da/dQei, as a function of —B-^- in Si. 



Then, the Green's function (G) becomes diagonal and is 
a solution of 



[E + h^/{2p)A - U{r)]GiE; r' , r) = -~5{r' - r). 



C/(r) = Us{r) + Uc{r), 



(A7) 
(A8) 



where fi is the reduced mass and Uc represents the 
Coulomb potential. Then, neglecting the off-diagonal 
couplings, the partial wave decomposition of S{E) can 
be expressed as 



JM Iy Jy riN.lN.jN 



(A9) 



where 



W{jN, Jy, J) = {2jN + l){iN\jQ\jY\f (AlO) 

for the closed shell, and 

'-'lYjY,nNlNjN y ' 

-\ j drdr'r^r'^[j*jM{p^,pK,eK,r)(f>l^i^j^{r) x 

(All) 

Here, (jy-nNiNON represents a radial wave function of a 
nucleon-hole state, and jjm is defined by the partial- 
wave decomposition of the distorted waves, written as 

x^'^*{PK,R)x^+\p.,R) = 

J2jJM{p.,PK,9K,r)Yfm. (A12) 

JM 
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G '=G'/ i „ , a , is a solution of the following equation 



The partial wave decomposition of the Green's function, 

,(f 2d ^(^+ 1), 

' I ^ nj — u yr 

\sir'-r), (A13) 



^2/i ^'dr"^ ' r dr 



-)+E^Uir)]G'^{E;r',r) 
1 



where G'^ is generally expressed as 

G' = Ay,irMr>)- 



(A14) 



Here, yi and j/2 are the regular and outgoing solu- 
tions of the differential equation, the left hand side of 
Ea. rn3l =0. and r< (r>) means the smaller (larger) one 



of r and r' . Factor A is normalized so as to satisfy 
Aiy2y{ - yiy'2)=i2t^)/{hr)\ 

Distorted waves of the pion and kaon were calculated 
using the eikonal appr oximat ion , assuming mb 
and (7i^+N~l4: mb [43 ■ The distortion reduces the mag- 
nitude of the spectrum. As described in Section V, the 
measured distortion seems to be stronger than that ex- 
pected from the eikonal approximation. The absolute 
strength of the distortion is still unclear. The strength 
of the distortion does not affect the spectrum shape very 
much. Thus, the magnitude of the spectrum is arbitrarily 
adjusted with a free parameter in the present spectrum 
shape analysis. 
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